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Abstract 

In this work we extend the range of applicability of a method re- 
cently introduced where coupled first-order nonlinear equations can be 
put into a linear form, and consequently be solved completely. Some 
general consequences of the present extension are then commented. 
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Despite of being less common, the systems governed by nonlinear equa- 
tions are of great importance and, nowadays, there is a growing interest in 
dealing with such systems pQ - [T3] . Unfortunately, as a consequence of the 
nonlinearity, in general we lose the capability of getting the complete so- 
lutions. Here we extend a demonstration [20] that is profitable for those 
systems in 1+1 dimensions, whose the second-order differential equations 
can be reduced to the solution of corresponding first-order equations, the so 
called Bolgomol 'nyi-Prasad-Sommerfield (BPS) topological solitons ^H]- I n 
this approach one can obtain a differential equation relating the two coupled 
fields which, once solved, leads to the general orbit connecting the vacua of 
the model. In fact, the "trial and error " methods historically arose as a con- 
sequence of the intrinsic difficulty of getting general methods of solution for 
nonlinear differential equations. About two decades ago, Rajaraman in- 
troduced an approach of this nature for the treatment of coupled relativistic 
scalar field theories in 1+1 dimensions. His procedure was model indepen- 
dent and could be used for the search of solutions in arbitrary coupled scalar 
models in 1+1 dimensions. However, the method is limited in terms of the 
generality of the solutions obtained and is convenient and profitable only for 
some particular, but important, cases JH]- Some years later, Bazeia and 
collaborators [TH] applied the approach developed by Rajaraman to special 
cases where the solution of the nonlinear second-order differential equations 
are equivalent to the solution of corresponding first-order nonlinear coupled 
differential equations. By the way, Bazeia and collaborators wisely applied 
their solution to a variety of natural systems, since polymers up to domain 
walls. Last year [20], one of us showed that many of the systems studied 
in can be mapped into a first-order linear differential equation and, 

as a consequence, can be solved in order to get the general solution of the 
system. In this work we are going to show that one can extend the range of 
applicability of that procedure, by taking into account situations where the 
potential of the system can not be written in terms of a superpotential as 
it is usually supposed to be necessary [Hj-j^HI- After that, we trace some 
comments about the consequences, and possible applications, coming from 
these general solutions. 

In order to deal with the problem, following the usual procedure to get 
BPS [TE] solutions for nonlinear systems, one usually particularizes [T§]-[2l)] 
the form of the Lagrangian density 

L = \{d,4>f + \{d,xf-V{^ X ), (1) 
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by imposing that the potential must be written in terms of a superpotential 
W ((/), x) such that 



(2) 



Then, the energy functional of the static configurations can be calculated 
straightforwardly, giving 



1 f°° 
En = - dx 



'a ,+ (s) ,+w ! +w ! - 



(3) 



which leads to 



E 



B 



dx 



(4) 



and finally, considering the minimal energy field configuration, one gets 



E B = \W(^ lXj )-W^ l ,Xi)l (5) 

where 0, and Xi are t ne are the % — th vacuum state of the model 22 = . 

In this case, one can easily see that solutions with minimal energy of 
the second-order differential equation for the static solutions in 1+1 dimen- 
sions, can be solved through the corresponding first-order coupled nonlinear 
equations 

^ = w^ tX )^ = w x (M, (6) 

where = ^ and W x = It is interesting to remark that the BPS so- 
lutions settle into vacuum states asymptotically. In other words, the vaccum 
states act as implicit boundary conditions of the BPS equations. 

In this work we are going to expand the class of nonlinear models with 
exact solutions by relaxing the condition over the potential to a less restrictive 
one. Namely we uses that 

V (0,X) = \{A{^x)f ] + \{B{^x)f ■ (7) 
As a consequence the energy functional can then be written as 
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E 



B 



dx 



dx I \ dx 



, dx dx , 



(8) 



Once again, it can be shown that the configuration which makes the 
energy minimal is that where the fields obey the first-order differential equa- 
tions 

f x=A(M ,f^ B{M . (9) 

This happens because, using the above equations we see that the energy 
functional becomes simply 



E 



B 



dx 



'dxY + fd£\ 2 ' 

dx I \dx I 



(10) 



which is positive definite, so rendering the minimal energy for the fields 
configuration obeying Q. 

Furthermore, in order to guarantee that the above equations obey the 
correct second-order differential equations 



d 2 



dx 2 



V<t> (4>, x) , 



dx 2 



V x (0, x) , 



(11) 



one must impose the following restriction over the, otherwise arbitrary func- 
tions, 

A x (<p, X ) = B (f) (<f ) ,x)- (12) 

Apart from this constraint and presenting a computationally more in- 
volved calculation of the topological energy, this more general situation can 
now be studied by using, for instance the usual trial-orbit approach [T§j-[25] 
or, by a direct solution of the orbit equation as done in j20]. In this work 
we choose the second line of thinking. The idea is to write the following 
equation 

(13) 



— = d = — 
A B ' 



where the spatial differential element is a kind of invariant. So, one obtains 
that 

# = A{<j),x) (u) 

d X B(<f>,x)' { ] 
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This last equation is, in general, a nonlinear differential equation relating 
the scalar fields of the model. Now, if one is able to solve it completely, 
the function <fi (x) can be used to eliminate one of the fields, so rendering 
the equations (jHJ) uncoupled. Finally, this uncoupled first-order nonlinear 
equation can be solved in general, even if numerically. 

At this point it is interesting to deal with a concrete example, in order to 
better expose the advantage of the present approach. Let us use, for instance 
the following functions A and B, as the potential defining ones, 



A{<j>,x) = 
B{<t>,x) = 



y 2 + y X 2 + 7i X + ^1 + ei X + Ci 



TT^ 2 + TTX 2 + 72 <P X + fa 4> + e 2 X + (2- 



(15) 



By requiring that the condition (J 12)) over the A and B functions be kept, 
we get the constraints: 



0i = 72; 7i = «2; e x = 8 2 



(16) 



If one try to construct a superpotential from the above functions A and 
B, in the spirit which it is usually done, one concludes that it is not possible 
to do so unambiguously, predominantly due to the difference between the 
one-field dependent terms. For instance in the above illustrative example 
one would obtain the following two alternative superpotentials 



(17) 



Wi = J Adfr W 2 = J Ad X , 
leading respectively to the following first-order BPS equations 

dtf> A ((f), x) 



and 



dx [£?(0,x)-(fx 2 + e 2X + C 2 ) 



^(0,x)-(f0 2 + 5i0 + d 



(18) 



d X B(<P, X ) ' (19) 

So, as we can see from the above, one can at least extend the usual BPS 
procedure through introduction of arbitrary functions depending on only one 
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of the interacting fields. In order to have an exact explicit example, we further 
restrict our parameters choice by putting 7x = (3 2 = $i = e i = e 2 = C2 = 0, 
which leaves us with 



#_A(0,x)_f0 2 + fx 2 + C 



dx B (0, x) Pi<Px 



(20) 



which, after defining the field a = (j) + ^"Cij leads us to the following linear 
differential equation [20] . 

do «i „ ,„ . 

(21) 

and finally to the following exact orbit equation 

= c oX ,3l +77^ X , 22 

2 Pi — «! 

with Co being the arbitrary integration constant. In this case the system 
potential is given by 

V (*, X) = f 4 + ^ 2 + f X* + + + f , (23) 

which is to be compared here with the one discussed by Boya and Casahorran 
some years ago [IB], and can be written, identifying their fields and ours 
through a = (ft as p = x, as 

Vbc (0, X) = -/ - ^ 2 + \x A + {I ^X 2 + + ~, (24) 

from which we can map the constants of the two models in order to see, if 
possible, in what situation they can me mapped. From an inspection of the 
above potentials we conclude that a\ = (1 = —^75, Pi = ± v^A, and 
also the following restriction over the original parameters in Vbc should be 
taken into account, 

d = V\ (2 V\ ± l) , / = 2 A, (25) 



which implies into a constraint over the these two last parameters such that 

d=f±4- 
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At this point it is important to stress that a system like that studied in 
[THj can not have their solutions reached by the BPS usual procedure |TU] - 
[2*5] . On the other hand, in their work ^H] Boya and Casahorran have made a 
careful analysis of the soliton stability, showing under what circumnstances 
they are stable under linear perturbations, which now can be extended to 
our case, at least when the above appropriate constraint holds. Furthermore, 
now in a similar fashion to what has been done in the case of superpotential 
procedure to get topological configurations [20], we can show that there are 
further solutions to those obtained in jTH] , at least when the above mentioned 
restrictions are kept. Those mentioned solutions have the special property 
of being kink-like for both fields, in contrast with happens with the already 
known ones ^H]- An important remark at this point, is that here the solutions 
are by no means obtained through a trial and error approach, instead, they 
are exactly solved. 

Let us consider now a bit more general model, where 71 = fa = £l = e 2 = 
£2 = 0. By using the same field transformation than in the previous case we 
get the following linear differential equation 

da 
dx 

which, after solved, leaves us with the exact orbit equation, connecting the 
two fields, 

_ |" 2e^ + 2ai e x x + ai Pi X 2 
a[X) ~[ ai (2h- ai ) 

for a.\ 7^ 2 Pi, and 

°- (X) = A + c i ( £ i + A X) 2 + ^^ft^ ln (ex +Pix), (28) 
z Pi pi 

for ai = 2 Pi; Cq and c\ are arbitrary integration constants. From now on, we 
substitute these solutions in one of the equations @, and solve it, so obtain- 
ing a general solution for the system. In general it is not possible to solve x 
in terms of <j) from the above solutions, but the contrary is always granted. 
Here we will substitute 4> (x) m t ne equation for the field x- However, for the 
case ai = 2 Pi it is very hard to get an exact solution and we will, as usually 
done in this type of situation [20] to the deal with the case where a\^2 P\. 



(aia + pix 2 + 2 six) 
(PiX + ei) 



(26) 



+ c (Ax + e 1 )ft , (27) 
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By doing the described substitution, one gets 



dx 
dx 



± (A x + ei) 



1 



2 e? + 2 ai ei x + ai ft X 2 



+ c (ftx + ei)" 1 



^ 2Ci 



a-i 



(29) 

In general we can not have an explicit solution for the above equation, 
but one can verify numerically that the solutions are always of the same 
classes. Notwithstanding, some classes of solutions can be written in closed 
explicit forms. First of all, we should treat the system when cq = 0, because 
in this situation we can solve analytically the system for any value of the 
parameters, apart from the case where a% = 2 ft as stated in above. 

In order to get an explicit exact solution in this situation, we make a 
re-scaling in the spatial variable, such that x = ft ?/, and note that starting 
from the Ricatti equation 



— — = a \(p + — ■ - 6 



(30) 



and making the field transformation ip = y/Ax 2 + Bx + C, we finish with 



dx 



B 



Now, performing the following identifications, 



(31) 



.4 



(2 ft -ax) 



B = 



(2 ft 



C 



Oil) 



2 
a i 



(2ft-ai) 



-Ci 



(32) 



we are able to connect the equations and (|3*T|) . which allows us to obtain 
the solution we are seeking for, just by using the solution of the Ricatti 
equation (|3U|). and substituting it in 



X± = 



2A 



or 



X 



B±^JB 2 -4A(C - tp 2 

. (v 2 - C) 



, A^O, 



B 



A = 0, B ^ 0. 



(33) 



(34) 



In this second hypothesis, which is not allowed for the set of parameters of 
the model we are considering, but which could appear in another system, we 
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must take into account that the starting Ricatti equation should be written 



as 



d( P _i_o f.Ji n e ± B ^ 




and performing the transformation tp = ^JB \ + we finish with 



(35) 



«fr ± ^ + Sj^ X + a (36) 

Our final task is to obtain the solution of ()30|) . In fact this a very well 
known equation but, for the sake of completeness, we map it too into a 
linear equation from which the solution of our original coupled nonlinear 
first-order differential equations will come after a sequence of mathematical 
manipulations. For this we perform a final transformation <p = ±Vo + i, 
where b = — C (A ^ 0), leaving us with 

| = «(l±2^), (37) 

from whose solution we obtain 

<p(y)=±Vb -, (38) 

with dq being the arbitrary integration constant. In order to avoid a singu- 
larity in the solution, we must choose the negative sign in the above solution. 
After that, it is not hard to verify that the final function can be expressed 
as the one having the usual kink profile, 

ip (y) = Vbtanh aVb(y — yo) . (39) 

However, once the field <p appears quadratically in the equations (p^j) 
and (|34|). the solution for this field will never have a kink profile. It will 
always have a lump behavior. This feature of the cq = case was already 
observed in and here we confirm that this behavior is still present in 
this more general case. Notwithstanding, we remark that if e\ 7^ 0, there is 
a remaining possibility in this more general system, which were not present 
in the simpler case studied in [20 . For this case, we can look for solutions 
in a special situation where the x field equation can be written in a simpler 
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polynomial form. This happens provided that B = \J A C, and in this case 
one gets 



dx ( 
and finally 



where x' 
solution 



B 



X+|) ^A X 2 + B X + C = a (x + ^J (VAx + Vc), (40) 



2A 



dx 
dy 



aVA( X '*-b'), 



(41) 



x + HvS + t) and & ' 



I (V^ ~ 1 with the kink - like 



X 



VF tank a V 'Ah' (y — y ) 




(42) 



It is not difficult to conclude that this last solution will lead to what we 
call type B kink [20], where both interacting fields x an d 4> present a kink 
profile, and the previous solution corresponds to that usually obtained in the 
literature jin],[2l]-[2H], what we call type A kink and is characterized by a 
lump form for the x field and a kink form for the (f> field. 

As observed in [20], for the case with cq ^ we can, at least for certain 
relations between the potential parameters, obtain analytical exact solutions 
for the type B kink solution. This can be achieved through a special fixing of 
the arbitrary integration parameter cq. For instance, in the case of the model 
on the screen in this work, we could choose cti = f3i and, after rearranging 
the terms in the equation (J40)) . conclude that the condition required in order 
to get a polynomial at the right-hand side of the equation for the field Xi 
could be granted by imposing the following condition over cq, 



B + c (3 l = sjA (C + c ei) 



(43) 



A similar analysis could be done for the case where ot\ = 4 0i which, despite 
of being more involved, lead us to the same conclusion about the possibility 
of generating type B kink through a convenient choice of cq. 

I this work we showed how to extend an approach recently introduced 
20J, in order to enlarge the class of nonlinear coupled field models which can 
have their solitonic configurations exactly solved in 1 + 1 dimensions. This 
allowed us to, at last under certain special conditions, include models like 
those studied by Boya and Casahorran [TH] in this class. As a bonus, we got 
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a kink configuration which was not reported in that work. Furthermore, it 
is important to assert that the class of models which we were able to deal 
here could not be considered through the usual procedure, based on the idea 
of having a superpotential Finally we should comment that similar 

models were very recently considered in the literature, concerning the study 
of Lorentz breaking models [2o1 127j. 
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